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Abstract 

When an atom is strongly coupled to a cavity, the two systems can exchange a single photon 
through a coherent Rabi oscillation. This process enables precise quantum-state engineering and 
manipulation of atoms and photons in a cavity, which play a central role in quantum information 
and measurement. Recently, a new regime of cavity QED has been reached experimentally where 
the strength of the interaction between light and artificial atoms (qubits) becomes comparable to 
the atomic transition frequency or the resonance frequency of the cavity mode. Here we show that 
this regime can strongly modify the concept of vacuum Rabi oscillations, enabling multiphoton 
exchanges between the qubit and the resonator. We find that experimental state-of-the-art circuit- 
QED systems can undergo two- and t/iree-photon vacuum Rabi oscillations. These anomalous Rabi 
oscillations can be exploited for the realization of efficient Eock-state sources of light and complex 
entangled states of qubits. 
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I. INTRODUCTION 


Light-matter interaction in the strong-coupling regime is a coherent reversible process in 
which a photon is absorbed and re-emitted by an electronic transition at a rate equal to the 
coupling energy divided by the Planck constant [1, 2], Reaching the light-matter strong¬ 
coupling regime has been a major focus of research in atomic physics and quantum optics for 
several decades and has driven the held of cavity quantum electrodynamics (cQED) [3, 4], 
The strong-coupling regime has been observed, both in the time and frequency domain, in a 
variety of systems [5], when an electronic transition is resonantly coupled to a cavity (optical 
resonator) and the coupling rate exceeds the rates of relaxation and decoherence of both 
the electronic transition and the held. Cavity QED ehects with individual qubits have been 
intensively studied in solid state systems by replacing natural atoms with artihcial atoms, 
such as quantum dots [6, 7] and Josephson circuits [8-11]. The strong-coupling regime, when 
reached with a single qubit, enables a high degree of manipulation and control of quantum 
systems [12]. For example, by exploiting the strong-coupling regime of cavity QED, the 
preparation and measurement of arbitrary quantum states in a completely controlled and 
deterministic manner has been achieved [13, 14] and “Schrodinger’s cat” states of radiation 
have been prepared and reconstructed [15, 16]. Basic steps in quantum information process¬ 
ing, including the deterministic entanglement of atoms and the realization of quantum gates 
using atoms and photons as quantum bits have also been demonstrated [2, 17, 18]. 

Recently a new regime of cavity QED, where the coupling rate becomes an appreciable 
fraction of the unperturbed frequency of the bare systems, has been experimentally reached 
in a variety of solid state systems [19-27]. In this so called ultrastrong-coupling (USC) 
regime, the routinely-invoked rotating wave approximation (RWA) is no longer applicable, 
and the antiresonant terms in the interaction Hamiltonian signihcantly change the standard 
cavity-QED scenarios (see, e.g., Refs. [28-36]). Although counter-rotating terms in principle 
exist in any real light-matter interaction Hamiltonian, their effects become prominent only 
in the USC limit [37]. Usually, light-matter USC is reached by coupling the resonator with a 
large number of molecules or more generally electronic transitions. USC with a single qubit 
has been achieved only by using superconducting circuits based on Josephson junctions, 
which exhibit macroscopic quantum coherence and giant dipole moments as artihcial atoms 
[20, 21, 38]. 
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One interesting featnre of the USC regime is that the number of excitations in the cavity- 
emitter system is no longer conserved, even in the absence of drives and dissipation. Mea¬ 
surements on a superconducting circuit QED system in the USC regime have shown clear 
evidence of this feature [21], Specihcally, tuning the qubit transition frequency (by adjust¬ 
ing the external flux bias threading the qubit), and measuring the cavity transmission, an 
anticrossing arising from the coupling between states with a different number of excitations 
has been observed. In particular, the measurements evidence the coupling of \g, 0, 0,1) and 
|e, 1,0,0), where the kets indicate the states of the qubit and of the hrst three resonator 
modes. Very recently it has been shown [39] that, when the frequency of the cavity held 
is near one-third of the atomic transition frequency, there exists a resonant three-photon 
coupling via intermediate states connected by counter-rotating processes. 

The resonant quantum Rabi oscillations, occurring when the atom and the cavity mode 
can exchange one excitation quantum in a reversible way, play a key role in the manipulation 
of atomic and held states for quantum information processing [12]. Here, we show that a 
system consisting of a single qubit ultrastrongly coupled to a resonator can exhibit anomalous 
vacuum Rabi oscillations where two or three photons are jointly emitted by the qubit into 
the resonator and re-absorbed by the qubit in a reversible and coherent process. We focus 
on the case of a hux qubit coupled to a coplanar-waveguide resonator, a system where the 
USC regime with a single artihcial atom has been demonstrated [21]. We hnd that this 
ehect can be observed at coupling rates of the same order of those already reached in these 
systems [21, 38]. 


II. DISSIPATION AND PHOTODETECTION IN THE USC REGIME 

In order to demonstrate multiphoton quantum Rabi oscillations in USC cavity QED, we 
calculate the time evolution of the mean output photon huxes and higher-order normally- 
ordered photon correlations. 

It has been shown that, in the USC regime, the usual normally-ordered correlation func¬ 
tions fail to describe the output photon emission rate and photon statistics. Clear evidence 
of this is that the standard input-output relations predict, even for a vacuum input and the 
system in the ground state, a hnite output photon flux proportional to the average number 
of cavity photons [32, 40, 41], i.e., (^)jut(^)^out(^)) oc (a^(f)a(f)), where and A~^^^{t) 
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are the positive- and negative-frequency components of the output held operator, while a 
and are the destruction and creation operators for cavity photons. A solution to this 
problem has been proposed in Ref. [32] . Considering for the sake of simplicity a single-mode 
resonator, it is possible to derive the correct output photon emission rate and correlation 
functions by expressing the cavity electric-held operator X = a -|- aMn the atom-cavity 
dressed basis. Once the cavity electric-held operator has been expressed in the dressed ba¬ 
sis, it has to be decomposed in its positive- and negative-frequency components and 
X~ [32]. Expanding the X operator in terms of the energy eigenstates \j) (with huj the 
corresponding eigenvalues) of the system Hamiltonian H, one hnds the relations 

A'+ = Y. X,k\])W ; -A- = (.Y+)', (1) 

j,k>j 

where Xj^ = + a)\k) and the states are labeled such that Uk > ojj for k > j. The 

resulting positive frequency output operator can be expressed as 

hu.(*) = A-(t) - ^C„.Y+, (2) 

where k, is the loss rate of the resonator due to the coupling to the external in-out modes 
and Cq is a constant proportional to the zero-point huctuation amplitude of the resonator 
[42]. Two aspects of these results are noteworthy: hrst of all, we note that in the USC 
regime, one correctly obtains X+lO) = 0 for the system in its ground state |0) in contrast 
to a 10) 7^ 0. Moreover, we notice that the positive-frequency component of X is not simply 
proportional to the photon annihilation operator d. As a consequence, for arbitrary degrees 
of light-matter interaction, the output photon flux emitted by a resonator can be expressed as 
"hout = k{X~X~^). Similarly, the output delayed coincidence rate is proportional to the two- 
photon correlation function {X~(t)X~(t + t)X^( t + t)X^( t)) . In quantum optics, it is well 
known that the signal directly emitted from the qubit is proportional to (d+d-). In circuit 
QED systems, this emission can be detected by coupling the qubit to an additional microwave 
antenna [14]. Indeed, in the USC regime the qubit emission rate becomes proportional to 
the qubit mean excitation number where are the qubit positive and negative 

frequency operators, dehned as C~^ = k>j C~ = {C^)\ with Cjk = 

(j|(d_ +d+)|fc). 

In order to properly describe the system dynamics, including dissipation and decoher¬ 
ence effects, the coupling to the environment needs to be considered. We adopt the master- 
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equation approach. However, in the USC regime the description offered by the standard 
quantum-optical master equation breaks down [43]. Following Ref. [41, 44], we write the 
system operators in the system-bath interaction Hamiltonian in a basis formed by the eigen¬ 
states of H. We consider T = 0 temperature reservoirs (the generalization to T 7 ^ 0 
reservoirs is straightforward). By applying the standard Markov approximation and tracing 
out the reservoir degrees of freedom, we arrive at the master equation for the density matrix 
operator p{t), 

P = ^ [P{t) , H] + £dampP(t) + ^4>Pit) , (3) 

where £dampP(t) = + r^^)^[b')(^l]p(^) with V[d]p = |(26p6t - pO^O - O^Op), 

describes dissipation effects arising from the resonator and qubit reservoirs. These cause 
transitions between eigenstates at rates 



( 4 ) 

P" = l\{j\»Ak)\\ 

( 5 ) 


where k and 7 are decay rates, here assumed to be spectrally constant, induced by the 
resonator and qubit reservoirs. Pure dephasing effects affecting the qubit are described in 
Eq. (3) by the last term C^p{t) = V[Y,j ^j\j){j\]p{t), where = ^/y^{j\a^\j), and 7 ^ is 
the pure dephasing rate. Note that only the most relevant diagonal contributions have been 
included. 

III. RESULTS 

Here we study a flux qubit coupled to a coplanar resonator in the USC regime. For suitable 
junctions, the qubit potential landscape is a double-well potential, where the two minima 
correspond to states with clockwise and anticlockwise persistent currents ±Jp [10, 11]. When 
the flux offset = 4>ext —*^ 0/2 = 0 , where $ext is the external flux threading the qubit and 
$0 is the flux quantum, the lowest two energy states are separated by an energy gap A. In the 
qubit eigenbasis, the qubit Hamiltonian reads Hq = hujqaz/2, where huq = /SP + {2Ip5^xy 
is the qubit transition frequency, which can be adjusted by an external flux bias. We note 
that the two-level approximation is well justified because of the large anharmonicity of this 
superconducting artificial atom. The resonator modes are described as harmonic oscillators, 
Hm = where is the resonance frequency, m is the resonator-mode index and 
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oL (®m) is the bosonic creation (annihilation) operator for the m-th resonator mode. We 
will consider A/2 and A/4 resonators. Then, the quantum circuit can be described by the 
following extended Rabi Hamiltonian 

H = H^ + J2[Hm + hgmXm{cos9a^ + sm9a^)]. ( 6 ) 

m 

Here, = dm + denote Pauli operators, Qm is the coupling rate of the qubit to 

the m-th cavity mode and the flux dependence is encoded in cos 9 = A/wq. The operator ax 
is conveniently expressed as the sum of the qubit raising ((j+) and lowering (d-) operators, 
which in the Heisenberg picture and for Qm = ^ oscillate as exp {iuJat) (negative frequency) 
and exp (—iu^t) (positive frequency) respectively. Thus, in contrast to the Jaynes-Cummings 
(JC) model [45], the Hamiltonian in Eq. (6) explicitly contains counter-rotating terms of the 
form a+dl^, d_dm, and d^dm- Considering only one resonator mode and a flux offset 

(5<ha; = 0, the Hamiltonian in Eq. (6) reduces to the standard Rabi Hamiltonian. 

A. Two-photon quantum Rabi oscillations 

We hrst consider the case of a flux qubit coupled to a A/2 superconducting transmission¬ 
line resonator with resonance frequencies = m'Kc/L^ where L is the resonator length. 
We use the qubit parameters A/h = 2.25 GHz and 2Jp = 630 nA, as in Ref. [21]. We are 
interested in the situation where the qubit transition energy is approximately twice that 
of the fundamental resonator mode: Wq ~ 2u\. We consider the qubit to be positioned at 
the center of the resonator, so that it does not interact with the resonator mode m = 2 
(see Fig. la). The other resonator modes are much higher in energy, detuned with respect 
to the qubit transition frequency by an amount signihcantly larger than the coupling rate: 
(a;])j—cjq) ~ {m — 2)ul, providing only moderate energy shifts for the coupling rates gm/(^m ~ 
0.2 considered here. We will thus only take into account the interaction of the qubit with 
the fundamental resonator mode. We diagonalize numerically the Hamiltonian from Eq. (6) 
and indicate the resulting energy eigenvalues and eigenstates as huji and |i) with i = 0,1 ,..., 
choosing the labeling of the states such that ujk > ooj for k > j. 

Figure lb shows the frequency differences = Ui — uq for the lowest energy states as 
a function of the qubit transition frequency, which can be tuned by changing the external 
flux bias The red dashed curves correspond to calculations obtained neglecting all 
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FIG. 1. (Color online) (a) Sketch of the distribution of the first three resonator modes (m = 1, 2, 3) 
of a transmission-line A/2 resonator. The resonance frequency of the first mode (blue solid curve) 
is set to u;f/27r = 4 GHz. The qubit is positioned at the center of the resonator, so that it does not 
interact with the m = 2 (red dashed curve) and is off-resonance to the m = 3 (green dotted-dashed 
curve) mode. The bare resonance frequencies of the second and third modes are = 2(jo\ and 
= 3a;)^, respectively. Qubit parameters are A/h = 2.25 GHz and 2/p = 630 nA. (b) Frequency 
differences Uifi = coi — coq for the lowest energy dressed states as a function of the qubit transition 
frequency cjq (which can be tuned by changing the external flux bias (54>a;) for the JC model (red 
dashed curves) and the extended Rabi Hamiltonian (blue solid curves). We consider a normalized 
coupling rate gi/uj\ = 0.15 between the qubit and the resonator. In both cases the ground state 
level is not displayed, (c) Avoided level crossing (blue solid curves) resulting from the coupling 
between the states |e,0) and \g,2) due to the presence of counter-rotating terms in the system 
Hamiltonian. The energy splitting reaches its minimum at Wq/27r ~ 7.97 GHz ~ 2{u}1/2tt). The 
anticrossing is not present in the JC model (red dashed lines), since it arises from the coherent 
coupling between states with a different number of excitations. 
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the counter-rotating terms (JC model). We observe a spectrum with two large-splitting 
anticrossings around cjq ~ bj\ which appear both in the dashed and the continuous curves. 
In the JC picture, they correspond to the resonant coupling between states with the same 
number of excitations. The lowest energy avoided crossing results from the coherent coupling 
of the states |e, 0) and If/,!), where g{e) indicates the ground (excited) state of the qubit 
and the second entry in the kets represents the photon number. When the splitting reaches 
its minimum, the resulting system eigenstates are 

^(K0>±|J.1))- (7) 

The higher-energy large avoided crossing in the plot corresponds to the second rung of the 
JC ladder, arising from the coupling of |e, 1) and 1(7,2). Only small quantitative deviations 
between the eigenenergies in the JC and in the extended Rabi model can be observed. 

When the counter-rotating terms are taken into account, the states |f) are no longer 
eigenstates of the total number of the excitation operator Wxc = -|- d+d-. For example, 

the system ground state can be expressed as a superposition |0) = Xln 
of bare states also involving nonzero excitations. Of course, when the normalized coupling 
gi/u:\ <C 1, only the coefficient c^g is signihcantly different from zero. Moreover, for 0 = 0 
parity is conserved [46, 47] and only states with an even number of excitations contribute 
to |0). The non-conservation of the total excitation number also affects the excited dressed 
states |j) = Xln'^gn.15') As ^ consequence, the dressed states |1) and |2) at the 

minimum splitting do not correspond to the simple JC picture of Eq. (7). 

The continuous line levels in Fig. lb also display a smaller amplitude avoided crossing 
when (Uq ~ 2ul. Observing that just outside this avoided-crossing region one level remains 
flat as a function of the flux offset with energy u ~ 2u:\ while the other shows a linear 
behavior with cjq, the splitting originates from the hybridization of the states |e,0) and 
\g-,2). This avoided crossing behavior is better shown in Fig. Ic, and the resulting states 
are well approximated by the states :^(|e,0) ± \g^2)). This splitting is not present in the 
RWA, where the coherent coupling between states with a different number of excitations is 
not allowed, nor does it occur with the standard Rabi Hamiltonian (0 = 0). 

Following the procedure described in Ref. [39], such a two-photon coupling between the 
bare states |e, 0) and |(7,2) can be analytically described by an effective Hamiltonian (see 
Appendix). As displayed in Fig. 2, the coupling between |e,0) and \g^2) can only occur via 
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FIG. 2. (Color online) Coupling between the bare states |e, 0) and \g, 2) via the intermediate states 
\g, 1) and |e, 1). The red dashed arrows and the solid black arrows indicate, respectively, the two 
different processes describing the two-photon resonant coupling between |e,0) and \g,2). 


the intermediate states \g, 1) and |e, 1). Indeed, if the system is initially prepared in the state 
|e, 0 ), two different processes can occur: either (i) the counter-rotating term a[az enables a 
virtual transition |e, 0 ) —>■ |e, 1 ), and then the term d|(j_ leads to the hnal transition to the 
state 1 ( 7 , 2 ); or (ii) the term a\a_ enables the transition |e, 0 ) —)■ \g^ 1 ) which is followed by 
the virtual transition l^f, 1 ) —>■ \g,2) induced by the term a\az- 

In order to obtain an analytical description of the effective coupling, we hrst reduce the 
extended Rabi Hamiltonian to the truncated Hilbert space composed of the bare states 
|e, 0), \g, 1), |e, 1), and \g, 2). The matrix form of the reduced Hamiltonian becomes 


n 


( 

2 

gi cos 6 
gi sin 9 

V 0 


gi cos 9 gi sin 9 0 ^ 

— 0 —a/ 2 ( 7 i sin6* 

0 + V^gi cos 9 

— ■\/2gism9 y/2gicos9 — ) 


( 8 ) 


where the order of columns and rows is |e, 0 ), \g, 1 ), |e, 1 ), and \g, 2 ). 

Near the two-photon resonance when Wq ~ 2,u:\, the intermediate states \g, 1) and |e, 1) 
can be adiabatically eliminated (see Appendix A), leading to the effective Hamiltonian 
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( 9 ) 


Hes = ^ cos(26')^ |e, 0 )(e, 0 | 

+ {^^1 “ ^ ^ cos( 20 )^ \g, 2 ) {g, 2 

-nf/"^(|e,0)(^?,2| + |f7,2)(e,0|), 
which describes the effective two-photon coupling between |e, 0 ) and \g, 2 ), with an effective 
two-photon Rabi frequency 

fiir> ^ (10) 

(Uq 

A key theoretical issue of the USC regime is the distinction between bare (unobservable) 
excitations and physical particles that can be detected. For example, when the counter¬ 
rotating terms are relevant, the mean photon number in the system ground state is different 
from zero: (0|ala|0) 7 ^ 0. However, these photons are actually virtual since they do not 
correspond to real particles that can be detected in a photon-counting experiment. Accord¬ 
ing to this analysis, the presence of an n-photon contribution in a specific eigenstate of the 
system does not imply that the system can emit n photons when prepared in this state. In 
order to fully understand and characterize this avoided crossing not present in the RWA, a 
more quantitative analysis is required. In the following, we therefore calculate the output 
signals and correlations which can be measured in a photodetection experiment. 

In order to probe the anomalous avoided crossing shown in Figs, lb and Ic, we consider 
the case where the qubit is directly excited via a microwave antenna by an optical Gaussian 
pulse. The corresponding driving Hamiltonian is 

Ha = £{t) cos{ujt)ax , (11) 

where £{t) = Aexp [—{t — to)^/(2r^)]/(rA/^). Here, A and r are the amplitude and the 
standard deviation of the Gaussian pulse, respectively. We consider the zero-detuning case 
by choosing the flux offset ( 5 <F 3 ; corresponding to the qubit frequency a;q/27r ~ 7.97 GHz, 
where the splitting in Fig. lb is at its minimum. The central frequency of the pulse has 
been chosen to be in the middle of the two split transition energies: 00 = (cua^o + ^ 2 , 0 )/‘2- 

The cavity output photon flux and the photon flux emitted by the qubit directly coupled 
to a microwave antenna are proportional to and respectively. Figure 3a 

displays the dynamics of these two quantities after the arrival of a vr-like pulse exciting the 
qubit described by the Hamiltonian (11). Galculations in Fig. 3 have been carried out in the 
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FIG. 3. (Color online) (a) Temporal evolution of the cavity mean photon number {X~X~^) (blue 
solid curve) and the qubit mean excitation number {C~C~^) (red dashed curve) after the arrival 
of a TT-like Gaussian pulse exciting the qubit (the black vertical line shows the wavepacket peak 
arrival time). The amplitude and the central frequency of the pulse are A/uj\ = 8.7 x 10“^ and 
u = (wa^o + W2 ,o)/ 2, respectively. After the arrival of the pulse, the system undergoes vacuum 
Rabi oscillations showing the reversible excitation exchange of two photons between the qubit 
and the resonator, (b) Time evolution of the zero-delay two-photon correlation function 
(dashed black curve) together with the intracavity photon number {X~X^) (blue solid curve). At 
early times they almost coincide. This perfect two-photon correlation is a signature that photons 
are actually emitted in pairs. Resonator and qubit damping rates are k/uiI = 1.8 x 10“^ and 
y/cjJ = 1.8 X 10““^, respectively. 


absence of pure dephasing ( 7 ,^ = 0). Results for 7 ^ 7 ^ 0 are shown in Fig. 4. Vacuum Rabi 
oscillations showing the reversible excitation exchange between the qubit and the resonator 
are clearly visible in Fig. 3a. The pulse time-width is not much narrower than the Rabi 
period, so the qubit excitation is partially transferred to the cavity during the pulse arrival. 
Therefore, the hrst peak of the qubit mean excitation number (Fig. 3a) is slightly lower than 
the second one. We observe that the mean intracavity physical photon number at its hrst 
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FIG. 4. (Color online) Effects of strong pure dephasing (dashed lines) on the dynamics of the mean 
photon number (a), the qubit effective population (b), and the two-photon correlation function (c). 
Calculations have been performed with the same parameters as in Fig. 3 with the addition of a pure 
dephasing rate 7 ^ = 300 7 . Solid lines display numerical results obtained in the absence of pure 
dephasing ( 7 ^ = 0). It can be observed that the physics of multiphoton vacuum Rabi oscillations 
is not significantly altered by the effects of pure dephasing. 

maximum is very close to two. This is a first hint that when the qubit is in the ground 
state the resonator mode acquires two photons. However, the output measured signals are 
proportional and not equal to {X~X^), so from this kind of measurements it is not possible 
to certify that the qubit and the resonator are actually exchanging two quanta. 

Figure 3b displays the time evolution of the zero-delay two-photon correlation function 
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= {X~(t)X~(t)X~^(t)X~^(t)) (dashed black line) together with the intracavity pho¬ 
ton nnmber {X~X^) (blue continuous curve) for comparison. At early times they almost 
coincide. This is a signature of perfect two-photon correlation: the probability of the system 
to emit one photon is equal to the probability to emit a photon pair. During the system 
evolution, higher values of the local minima of {X~X^) are observed due to the decay from 
the two-photon state to the one-photon state (out of resonance with respect to the qubit) 
caused by the photon escape from the resonator. However, the two-photon correlation 
function {X~(t)X~(t)X~^{t)X~^(t)) goes almost to zero every time the qubit is maximally 
excited. This different behavior in Fig. 3b indicates that the qubit does not absorb single 
photons but only photon pairs. The period of a complete population oscillation is 2ti 
where is the minimum splitting in Fig. Ic. 

Ordinary quantum vacuum oscillations have already been demonstrated in circuit-QED 
systems {e.g., [10, 48]). The dynamics observed in Fig. 3 can also be obtained by hrst 
preparing the qubit in its excited state by employing a vr pulse. Then the vr pulse can be 
followed by a shift pulse which brings the qubit into resonance with the resonator for the 
desired duration in order to observe the coupled dynamics as in Ref. [48]. If the shift pulse 
has a duration = 2tt the Fock state n = 2 is directly generated. After the switch- 
off of the shift pulse, the qubit is out of resonance with the resonator and the Fock state can 
escape from the cavity through an input-output port and be detected. Hence two-photon 
Rabi oscillations can be exploited for fast and efficient generation of two-photon states. 

The influence of strong pure dephasing effects is shown in Fig. 4. Calculations have been 
performed with the same parameters used in Fig. 3 with the addition of a pure dephasing 
rate = dOOy. Figure 4 compares the dynamics of the mean photon number (a), the qubit 
effective population (b), and the two-photon correlation (c) in the absence (continuous 
curves) and in the presence (dashed curves) of pure dephasing. The hgure shows that 
strong pure dephasing does not signihcantly alter the physics of multiphoton vacuum Rabi 
oscillations and the main effect of dephasing is to make the excitation pulse less effective. 


B. Three-photon quantum Rabi oscillations 

Very recently, it has been shown [39] that the strong coupling of a single qubit with 
three photons can be achieved in the USC regime when the frequency of the cavity held is 
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near one-third of the atomic transition frequency. In this case, parity-symmetry breaking 
is not required and this effect can occur also at 0 = 0. Hence, it could be observed also in 
systems like natural atoms or molecules displaying parity symmetry. One possible problem 
with this configuration is that the qubit can also interact resonantly with the one-photon 
state of the m = 3 mode of the resonator. In this case the qubit would interact with both 
one- and three-photon states. We show that the undesired one-photon resonant coupling 
of the qubit with a resonator mode can be avoided by considering a A/4 resonator, whose 
resonance frequencies are = (2m — l)7rc/2L, with m = 1, 2, 3 ... (see Fig. 5a). As shown 
in Fig. 5a, the qubit is positioned so that it does not interact with the mode m = 2 with 
resonance frequency = 3a;);. The qubit parameters are A/h = 4.25 GHz and 2/p = 630 
nA. In the present case we are interested in the situation where the qubit transition energy 
is approximately three times that of the fundamental resonator mode: cuq ~ 3a;). The mode 
m = 3 has resonance frequency o;^ = 5a;), which is much larger than a;q. Hence also in 
this case, we can consider the interaction of the qubit with only the fundamental resonator 
mode. 

Figure 5b displays the frequency differences Uifl for the lowest energy states as a function 
of the qubit transition frequency. The red dashed curves corresponds to calculations obtained 
neglecting all the counter-rotating terms (JC model). We observe a spectrum with two large- 
splitting anticrossings which appear only in the continuous curves plus a smaller avoided 
crossing magnihed in Fig. 5c. The lowest energy splitting corresponds to a two-photon 
vacuum Rabi splitting. When it reaches its minimum (at a;q ~ 2a;)), the corresponding 
hybridized states are analogous to those whose dynamics has been described in Fig. 3. They 
can be approximately expressed as 

2j(|e,0)±|9.2». (12) 

The second avoided crossing at higher energy corresponds to the second rung of the two- 
photon Rabi ladder and the corresponding approximated hybridized states (at the minimum 
splitting) are 

/|(k,l)±l9.3». (13) 

The third smaller splitting, occurring at a;q ~ 3c;), corresponds to a three-photon vacuum 
Rabi splitting. Here, a single qubit is resonantly coupled with a three-photon state, resulting. 
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FIG. 5. (Color online) (a) Sketch of the distribution of the first three resonator modes (m = 1, 2, 3) 
of a transmission-line A/4 resonator. The resonance frequency of the first mode (blue solid curve) 
is wj = 4 GHz. The qubit is positioned so that it does not interact with the m = 2 mode (red 
dashed curve) with resonance frequency = 2>uj\. The mode m = 3 (green dotted-dashed curve) 
has resonance frequency = buj \, which is much larger than uiq so that only the interaction of the 
qubit with the fundamental resonator mode can be considered. Qubit parameters are A/h = 4.25 
GHz and 2Ip = 630 nA. (b) Frequency differences uji^ = iVi — uq for the lowest dressed energy 
states as a function of the qubit transition frequency cjq (which can be tuned by changing the 
external flux bias (5<I>2;) for the JC model (red dashed curves) and the extended Rabi Hamiltonian 
(blue solid curves) explicitly containing counter-rotating terms. We consider a normalized coupling 
rate gi/u^l = 0.25 between the qubit and the resonator. The spectrum shows two large-splitting 
anticrossings, which appear only in the continuous curves, plus a smaller avoided crossing which 
is magnified in (c). (c) Three-photon vacuum Rabi splitting (blue solid curves) resulting from the 
coupling between the states |e,0) and Iff, 3) due to the presence of counter-rotating terms in the 
system Hamiltonian. The energy splitting reaches its minimum at ujq/2'K ~ 11.89 GHz ~ 3(w^/27r). 
The anticrossing is not present in the JC model (red dashed curves), since it arises from the coherent 
coupling between states with a different numbq^^of excitations. 
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FIG. 6. (Color online) (a) Temporal evolution of the cavity mean photon number {X~X'^) (blue 
solid curve) and the qubit mean excitation number {C~C~^) (red dashed curve) after the arrival of a 
TT-like Gaussian pulse exciting the qubit (the black vertical curve shows the wavepacket center time). 
The amplitude and the central frequency of the pulse are Ajuf-y = 9.4 x 10“^ and w = (a;4_o+^^3,o)/2, 
respectively. After the arrival of the pulse, the system undergoes vacuum Rabi oscillations showing 
the reversible excitation exchange between the qubit and the resonator. The fact that the mean 
intracavity physical photon number at its first maximum is very close to three is a first signature 
that, when the qubit is in its ground state, the resonator mode is in a three-photon state, (b) Time 
evolution of the zero-delay three-photon function (dashed green curve together with the 

intracavity photon number {X~X^) (solid blue curve). The first peak value of the three-photon 
correlation function is approximately two times higher than that of the mean photon number, a 
signature of an almost-perfect three-photon correlation. Parameters for resonator and qubit losses 
are the same as in Fig. 3. 
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Figure 6 displays the system dynamics after the arrival of a vr-like pulse exciting the 
qubit described by the Hamiltonian (11). Specihcally, Fig. 6 a shows the time evolution of 
{X~X^) and Calculations have been carried out in the absence of pure dephasing 

(70 = 0). Vacuum Rabi oscillations showing the reversible excitation exchange between the 
qubit and the resonator are clearly visible. We observe that the mean intracavity physical 
photon number at its first maximum is very close to three. This is a first hint that when 
the qubit is in the ground state the resonator mode is in a three-photon state. The period 
of a complete population oscillation is 27 ^where 211 is the minimum splitting in 
Fig. 5c. 

Figure 6 b displays the time evolution of the zero-delay three-photon correlation func¬ 
tion = {X~{t)X~{t)X~{t)X~^{t)X~^{t)X~^(t)) together with the intracavity photon 

number {X~X~^) for comparison. At early times the peak values of are approxi¬ 

mately two times higher than those of the mean photon number {X~{t)X~^{t)). This is a 
specihc feature of three-photon Fock states and indicates an almost perfect three-photon 
correlation. We observe that G^^\t) at early times reach a peak value slightly beyond 6 . 
This indicates that the system has a nonnegligible probability to emit more than three 
photons. This is conhrmed by the presence of a non-zero four-photon correlation function. 
Analyzing the different transitions contributing to G^^\t), we can attribute this effect to 
additional low-frequency transition |4) —)■ |3). These transitions between Rabi-split states 
occurs when parity symmetry is broken, [49] and in this case produces a 4-photon cascade: 
|4) —)• 13) —?• 12) —)■ |1) —)■ |0). This small contribution cannot be observed if its low fre¬ 
quency is outside the frequency-detection window. Analogously to the two-photon case (see 
Fig. 4), pure dephasing does not signihcantly affect the dynamics of multiphoton quantum 
Rabi oscillations (plot not shown). 


C. Generation of entangled GHZ states 

Standard vacuum Rabi oscillations have been exploited for the realization of atom-atom 
entanglement. Here we show that the multiphoton Rabi oscillations can be directly applied 
to the deterministic realization of more complex entangled states. As a hrst application 
we discuss the deterministic realization of multi-atom Greenberger-Horne-Zeilinger (GHZ) 
states [50] by using only one resonator. The GHZ states lead to striking violations of local 
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realism and are an important resource for quantum information processing [51], quantum 
cryptography [52] and error correction protocols [53]. Superconducting circuits have been 
used to study GHZ states (see, e.g., Ref. [54, 55]). 

Consider a resonator coupled to qubit 1 in the USC regime where two-photon vacuum 
Rabi oscillations can occur. The resonator also interacts in the strong (not ultrastrong) 
coupling regime with two additional qubits (2 and 3). Although the coupling rates between 
the resonator and the qubits are fixed, the qubit-resonator interaction can be switched 
on and off by adjusting the qubit frequencies [14]. The protocol is simple and consists 
of three steps, one for each qubit. We start by exciting the ultrastrongly-coupled qubit 
with a TT-pulse. Then, by changing the flux offset (at time t = 0), we drive it into reso¬ 
nance with the two-photon state of the resonator (Fig. Ic). The system state at time t is 
\^jJ) = cos (r2gg^^^t)|e, g, g, 0) -l-sin(r 2 gg^'^^f)| 5 f, g, g, 2). We let the qubit interact for a 7 r /2 Rabi 
rotation so that the resulting state is (|e, 5 ^,( 7 , 0) -|- \g, g, g,2))/\/2. We then drive the qubit 
back out of resonance, stopping the Rabi rotation. The second step consists of driving qubit 
2 into resonance with the 1 -photon state of the hrst resonator mode for a tt rotation time so 
that the resulting state is: (|e, g, g, 0) — \g, e, g,l))/ a/2. For the third step, we similarly drive 
the third qubit into resonance with the hrst resonator mode for a vr rotation. The resulting 
state is: {\e,g,g,0) + |( 7 , e, e, 0))/a/2. At this point, the photon state can be factored out, 
leaving us with a three-qubit GHZ-like entangled state. A more conventional GHZ state can 
be obtained by sending a further vr pulse to the hrst qubit, so that the resulting hnal state is 
{\g, g, g, 0) -|- |e, e, e, 0))/ a/2. This procedure can be easily generalized to four qubits or more. 
In general, if n-photon Rabi oscillations are achieved, then n -|- 1 qubit GHZ states can be 
produced. We notice that this protocol does not need the initial synthesis of photonic or 
atomic superposition states [56, 57]. 


IV. CONCLUSION 

We have investigated vacuum Rabi oscillations in the USC regime. According to the 
Jaynes-Cummings model, the qubit and the resonator can exchange a single excitation 
quantum through a coherent Rabi oscillation process. Such Rabi oscillations play a key 
role in the manipulation of atomic and held states for quantum information processing [ 12 ]. 
Our theoretical predictions show clear evidence for physics beyond the Jaynes-Cummings 
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model and extend the concept of quantum Rabi oscillations. We find that multiphoton 
reversible exchanges between an individual qubit and a resonator can be observed in the 
use regime. Specihcally, we have shown that experimental state-of-the-art circuit-QED 
systems can undergo two- and three-photon vacuum Rabi oscillations. Still increasing the 
coupling rate, a higher number of photons can be exchanged with the qubit during a single 
Rabi oscillation. These anomalous Rabi oscillations can be exploited for the realization of 
efficient Fock-state sources of light, and for the implementation of novel protocols for the 
control and manipulation of atomic and held states. 
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Appendix A: Analytical derivation of the two photon-qubit effective Hamiltonian 


In this Appendix, we derive the analytical expression for the effective Hamiltonian in 
Eq. (9), describing the two-photon coupling between the states |e,0) and 1(7,2). We start 
from the reduced Hamiltonian in Eq. (8) and then move to the rotating frame with frequency 
Ci;q/2, obtaining the transformed reduced Hamiltonian 



^ 0 

gi sin 9 

gi cos 9 

0 ^ 


gx sin Q 

ul 

0 

\/2gi cos 9 

h 

gi cos 9 

0 

U\ - Wq 

— \/2gi sin6' 


\ 0 

V^gi cos 9 

— \/2gi sin^ 

2ul - Wq . y 


(Al) 


Now, the order of columns and rows is |e, 0), |e, 1), \g,l), and \g,2). After the transfor¬ 
mation, an arbitrary state of the system in this truncated Hilbert space can be denoted as 
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(cl, c2, c3, c4)^ and the Schrodinger equation with Hamiltonian H'^ gives 


ici = {gi sin 0)02 + {gi cos 0)03 (A 2 ) 

ic2 = oj\c2 + {gi sin 6 )ci + (cos 6*)c4 (A3) 

ics = (!jj\ — a;q)c 3 + {gi cos6*)ci — (v^(7i sin6*)c4 (A4) 

iCi = {2oj[ — a;q)c4 + {\/2gi cos 6)02 — {y/2gi sin 6)03 . (A5) 


For gi/ul 1, the adiabatic elimination in Eqs. (A3) and (A4) can be applied [39] and the 
coefficients C 2 and C 3 can be approximated as 


C 2 ~ — — (sin 6 *ci + \/2 cos 6 * 04 ) 




C3 


9i 


-{cosOci — a /2 sin 6 ^ 04 ). 


(A 6 ) 

(A7) 


- Wq) 

The coupled equations for ci and C 4 are obtained substituting these results in Eqs. (A2) and 
(AS): 

^ ^ g;(a.,W( 2 a)-^;) x/ 2 g;^.in( 2 a) 

a;K^-a;q) 2ul{ul - u^) ^ ^ 

-2^?(a;^-a;qcos2 0)+a;j [2(^)2-3^a;q + a;2] ^^2 m( 20 ) 

^C 4 ^ ;-r-^C 4 H- ^ -—Cl . A9 

^(a;];-Wq) 2u\{u\-u:^ 

Considering the near-resonant case u}\ ~ C(;q/2, transforming back to the laboratory frame 
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FIG. 7. (Color online) Comparison between the minimum energy splitting 2Hgg^^^/a;q obtained 
analytically (red solid line) and numerically (blue points) as a function of gi/ojq for 0 = 7r/4. 


and keeping only the g\ dependence terms in the diagonal elements, the effective Hamiltonian 
in Eq. (9) is obtained. According to the effective Hamiltonian (9), the ratio of the minimum 
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splitting at the avoided crossing (see Fig. Ic) to the qubit frequency Wq is given by 

>(2ph) / _ \ 2 


^^^=4y2sin(20)(^) . 

Wq V^q/ 


A comparison between analytical and numerical results for the minimum energy splitting 


2Q 


(2ph) 

eff 


/uq is shown in Fig. 7 as a function of gi/ojq 
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